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QCD Thermodynamics on lattice 

p. Petreczky ^ 

'^Physics Department and RIKEN-BNL, Brookhaven National Laboratory, Upton, NY 11973 USA 

I discuss recent developments in lattice QCD thermodynamics on the nature of the transition at finite tempera- 
ture and density, equation of state, screening of static charges and meson spectral functions at high temperatures. 



1. Introduction 

It is expected that strongly interacting matter 
shows qualitatively new behavior at temperatures 
and/or densities which are comparable or larger 
than the typical hadronic scale. It has been ar- 
gued by Hagedorn that production of hadronic 
resonances may lead to a limiting temperature 
above which hadronic matter can no longer exist 
The corresponding temperature can be esti- 
mated to be 174MeV Wr On the basis of asymp- 
totic freedom one would expect that dominant de- 
grees of freedom at very high temperature and/or 
densities are quarks and gluons which are no 
longer subject to confinement 3 . Such qualita- 
tively different behavior of usual hadronic matter 
with confinement and chiral symmetry breaking 
and this new state of matter (usually called Quark 
Gluon Plasma) would suggest a phase transition 
at some temperature. The existence of such phase 
transition was first shown in the strong coupling 
limit of QCD 0], followed by numerical Monte- 
Carlo studies of lattice SU(2) gauge theory which 
confirmed it 

Since these pioneering studies QCD at finite 
temperature became quite a large subfield of lat- 
tice QCD (for recent reviews on the subject see 
Refs. |t)l7l8j ). One of the obvious reasons for this 
is that phase transitions can be studied only non- 
perturbatively. But even at high temperatures 
the physics is non-perturbative beyond the length 
scales 1 / (g^ T) {g^ (T) being the gauge coupfing) 
0. Therefore lattice QCD remains the only tool 
for theoretical understanding of the properties of 
strongly interacting matter under extreme con- 
dition which is important for the physics of the 
early universe and heavy ion collisions. In this re- 
view I will discuss recent developments on studies 
of the QCD phase transition at finite temperature 
and density, equation of state, the problem of in- 
medium modifications of inter-quark forces and 



meson spectral functions. 

2. QCD transition at finite temperature 

One of the basic questions we are interested 
in is what is the nature of the transition to the 
new phase or if there is one at all, what is the 
temperature where it happens ^ and what is the 
equation of state. In addition one may wonder 
what drives the transition in QCD, which I will 
discuss in section 5. 

In the case of QCD without dynamical quarks, 
i.e. SU(3) gauge theory most of these ques- 
tions have been answered. It is well established 
that the phase transition is 1st order JU]. Us- 
ing standard and improved actions the corre- 
sponding transition temperature was estimated 
to be Tc/Vct = 0.632(2) or 265{l)MeV (as- 
suming a — (420Afey)^ for the string tension). 
RG-improved actions gave similar estimates for 
Tcl\f^ CU- When is expressed in units of 
the Sommer scale tq comparison between differ- 
ent action can be done in a quite transparent way 
and gives r^T^ = 0.7498(50) [E]. The equation 
of state of SU(3) gauge theory were also extrap- 
olated to the continuum limit |13I14| . 

The situation for QCD with dynamical quarks 
is much more difficult. Not only because the in- 
clusion of dynamical quarks increases the compu- 
tational costs by at least two orders of magnitude 
but also because the transition is very sensitive to 
the quark masses. Conventional lattice fermion 
formulations break global symmetries of contin- 
uum QCD (e.g. staggered fermion violate the fla- 
vor symmetry) which also introduces additional 
complications. 

Since simulations at the physical value of the 
quark masses are extremely expensive one usu- 



will talk about QCD finite temperature transition irre- 
spective whether it is a true phase transition or a crossover 
and Tc will always refer to the corresponding temperature. 
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Figure 1. The QCD phase diagram on the plane 
of iTiu^d and nis. 

ally studies the QCD phase diagram on the mass 
plane of u,d and s-quarks. This is schematically 
shown in Fig. ^ As mentioned before for in- 
finitely heavy quarks there is a deconfinement 
transition which is well known to be 1st order. 
The 1st order nature of the transition persists 
for finite values of the quark masses and termi- 
nates with a critical line which belongs to Z{2) 
universality class at quark masses corresponding 
to pseudo-scalar meson mass of about 2.5GeV 
|15|. For the case of three fiavors of light quarks 
universality arguments suggest a 1st order chi- 
ral transition which was indeed observed in lat- 
tice simulations with both standard and improved 
(p4) staggered fermion actions |16ll7j . The re- 
gion of the 1st order transition ends at some crit- 
ical line which belongs to Z{2) universality class 
|lf)ll7| . Note, however, that standard and im- 
proved fermions give different estimates for the 
value of quark (pion) mass where the 1st order 
region ends, ~ 270MeV for standard action 
Ua and = 67(18)MeF for p4 action 
The transition temperature in the chiral limit for 
iVj = 4 was estimated to be = 154(8) [T^ . 

For two massless flavors one expects a 2nd or- 
der transition with 0(4) universality class. Sim- 
ulations with staggered fermions (both standard 
and improved) did not confirm this expectations. 
On contrary a recent analysis with standard stag- 
gered action on large lattices suggests that tran- 
sition maybe weakly 1st order 19 . On the other 
hand simulations with Wilson fermions show fi- 



nite size scaling consistent with 0(4) exponents 
|2()j . Despite these disagreements the chiral crit- 
ical temperature for two fiavor QCD turns out 
to be consistent for Wilson and staggered for- 
mulations: for Wilson fermions one has Tc = 
171(4)Mey 1201 1 while for improved staggered 
(p4) Tc = 173{7)MeV W ( both refer to TV* = 4 
lattice). First attempts to study the finite tem- 
perature transition in two flavor QCD with chiral 
fermions which have all the symmetries of contin- 
uum QCD have been done in Ref. |21| . 

The most interesting question is about the na- 
ture of the transition in the physical case of two 
light and one heavier quarks and where the phys- 
ical point in Fig. ^ is located. Simulations 
for this case have been performed using stan- 
dard |,22.23j and improved staggered formulations 
|24I25I26I27| . The most recent calculations with 
standard staggered fermions were performed on 
Nt = 4: lattices and quark masses very close to 
physical ones corresponding to the lightest pion 
mass of about WOMeV ( m^/nip = 0.188(2)). 
Calculations with improved staggered (Asqtad) 
action were done using three different lattice spac- 
ings (Nt) and several quark masses, with light- 
est quark mass being only 2.5 times heavier than 
the physical value. Both studies have found no 
evidence for a true phase transition but only a 
crossover. The crossover temperature was esti- 
mated to be Tc = 164(2)Mey (the errors is only 
statistical and tq = 0.5 fm was assumed) for stan- 
dard staggered action and 169(12)(4)MeV 
for Asqtad The later is an extrapolated 

value. The continuum and chiral extrapolation 
was done with Tc{Nt) = Tc+ci{m.„ /mpY+C2/Nf 
Ansatz with d — 1.08 motivated by 0(4) scal- 
ing. The coefficient C2 turns out to be quite small 
as expected; in Asqtad action the anomalously 
large agO?' corrections are removed. The second 
error is the estimate of Tc is the uncertainty in the 
scale (lattice spacing) l^^. If expressed in units 
of the Sommer scale ro both standard action and 
Asqtad action give a consistent estimate for Tc- 
When using the most recent value tq = 0.467/m 
[2H| and rescaling the results one gets the value 
Tc = n&{2)MeV for the standard action. 

It seems that different fermion actions give very 
similar estimate for the transition temperature. 
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Figure 2. The chiral susceptibility in three flavor QCD for A^t = 4, 6 and 8. 



which interestingly enough appears to be very 
close to the old prediction from the resonance 
gas model. Why are there so large discrepan- 
cies concerning the order of the phase transition? 
As mentioned before staggered fermions violate 
the flavor symmetry. This results in the fact that 
there is only one Goldstone pion, the other pions 
remain massive in the chiral limit for any finite 
lattice spacing. For typical lattice spacings used 
in thermodynamic studies (a — 0.3 — 0.2 fm) the 
splitting between Goldstone and non-Goldstone 
pions can be quite large, so the lightest non- 
Goldstone pion turns out to be heavier than the 
kaon. Clearly it is important to improve the flavor 
symmetry of staggered fermions when one wants 
to make statement about the nature of the tran- 
sition at the physical point. This is usually done 
by introducing fat links in the fermion hopping 
term (see e.g. Refs. |18I24| ') . Typically the fat 
links which are used are the sum of a single link 
and staples and are not elements of the SU(3) 
group. In staggered HYP action the fat links are 
elements of SU(3) gauge group ^Hl- It has been 
shown that projecting the fat links onto SU(3) 
group substantially improves the flavor symme- 
try (29,. Finite temperature simulations with 
HYP actions gave a number of interesting re- 
sults . For QCD with four light flavors us- 
ing standard staggered action the transition was 
found to be first order in accord with universal- 
ity arguments '3 0I31| . The only problem is that 
the critical quark mass where the transition turns 
into crossover was found to be strongly dependent 
on the lattice spacing it decreases by more 



than a factor of three when going from A^t = 4 to 
Nt = 6. In Ref. 32 it was found that introducing 
unprojected fat links make the transition weaker, 
while with the use of HYP smeared fat links the 
first order transition is washed out completely. In 
fact for HYP action, no first order transition was 
observed even for very small quark masses both 
for Nt — 4: and Nt — 6 lattices 32 . There are 
also simulations with HYP action at = 4 in 2 
flavor as well as at TVj = 4 and 6 in 2 -f 1 flavor 
QCD which indicate that the crossover is much 
smoother than for standard staggered action for 
the same quark mass and same lattice spacing 
j33l34| . Can improvement of flavor symmetry in 
the staggered fermion action smoothen the tran- 
sition and eventually wash out any phase transi- 
tion for the quark masses where standard action 
predict a phase transition? In this respect it is 
interesting to study the smoothness (sharpness) 
of the crossover as function of the lattice spac- 
ing. Typically susceptibilities of order parameter 
like quantities, e.g. chiral condensate ipip, show a 
pronounced peak for a rapid crossover. For a true 
phase transition the height of the peak diverges 
in the thermodynamic limit, while for a crossover 
it approaches a finite value. The MILC collabora- 
tion has calculated the susceptibility of the chiral 
condensate for three flavor QCD at three differ- 
ent lattice spacings 24 25 26 27 . The results are 
shown in Fig. |21 One can see a pronounced peak 
for Nt — 4, however, the peak height relative to 
the low temperature value of the chiral suscepti- 
bility decreases with increasing Nt- This means 
that as the flavor symmetry gets restored by go- 
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ing to finer lattices (smaller Nf) the crossover be- 
comes smoother. Note that the chiral susceptibil- 
ity in Fig. 13 is given in lattice units. To convert 
to physical units, say units of the data should 
be multiplied by which will make the chiral 
susceptibility look very different for different Nf. 
Thus it is likely that the first order transition ob- 
served so far with staggered action is just a lattice 
artifact. 

3. Phase diagram at /i > 

For a long time there was no considerable 
progress in studying the phase diagram of QCD 
at finite chemical potential. This is because the 
complex fermion determinant makes importance 
sampling based simulations impossible. The 
multi-parameter re-weighting method by Fodor 
and Katz was a partial breakthrough which al- 
lowed to study the phase diagram of QCD also at 
finite chemical potential |35I36) . Soon after also 
other methods based on Taylor expansion around 
/I = |;-{7l88j and analytic continuation to imagi- 
nary chemical potential |39l4()l4lj were suggested. 
All three methods were discussed in detail by S. 
Katz at Lattice 2003. Therefore I will emphasize 
mostly the most recent results leaving technical 
details aside. 

In the previous section we have seen that all 
available lattice data suggest that there is no 
phase transition in real QCD but only a smooth 
crossover. In this case a very interesting phase di- 
agram in the plane of temperature and chemical 
potential (T, /i) has been conjectured. At large 
chemical potential a 1st order transition is ex- 
pected and thus if there is a crossover at fi = the 
phase transition line should end at some critical 
point Te, ^J■E■ The important question is whether 
we can locate this end-point on the lattice. The 
first attempt has been made by Fodor and Katz 
using multi-parameter re- weighting [Hfij . At each 
value of /i the pseudo-critical point in the gauge 
coupling (3 was located using Lee- Yang zeros, ze- 
ros of the partition function for complex /3 which 
are closest to the real axis. In the crossover re- 
gion Lee- Yang zeros have non-vanishing imagi- 
nary parts in the thermodynamic limit. For a true 
phase transition, on the other hand, the imagi- 
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Figure 3. The phase diagram at finite chemi- 
cal potential obtained using multi-parameter re- 
weighting 23 . 



nary part of Lee- Yang zeroes vanishes in the in- 
finite volume limit. Thus inspecting the finite 
size scaling of the imaginary part of Lee Yang ze- 
ros one can determine whether there is a phase 
transition and eventually determine its order; the 
real part of Lee- Yang zero gives the location of 
the critical /3 (or equivalently the temperature). 
Using this procedure with standard staggered ac- 
tion, quark masses about three times larger than 
the physical one and lattices 4^, 6^x4, 8'^ x 4 they 
estimated the location of the critical end-point at 
^ 172{3)MeV, ^1 = 725{35)MeV ^ (note 
that the baryon chemical potential fiB = -^fJ- with 
fi being the quark chemical potential). Very re- 
cently they extended their studies by including 
larger lattice volumes and values of the quark 
masses very close to the physical ones The 
result of their new analysis is summarized in Fig. 
13 The endpoint was found to be located lo- 
cated at = 162(2)Mey, = 360(40)MeV. 
It is evident that the end-point is very sensitive 
to quark mass. Calculations using the method 
of analytic continuation in 2-1-1 flavor QCD ,3^ 
with the same quark masses as used in Ref. |3(i| 
give results which are in very good agreement 
with findings of Fodor and Katz for the position 
of the pseudo-critical line. There are also cal- 
culations based on Taylor expansion in 2 flavor 
QCD with somewhat heavier quark masses which 
gives qualitatively similar results. In the case of 
3 flavor QCD using re-weighting with Taylor ex- 
pansion the critical endpoint was estimated to be 



fi% = 120(30)Mey for 



UOMeV [ini. 
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4. Equation of state 

Lattice calculations of the equation of state are 
difficult for at least two reasons. First, the calcu- 
lation of the energy density and pressure needs 
subtraction of the zero temperature (vacuum) 
part. This makes computational expenses grow 
with lattice spacing as a^^^ compared to 
growth for other quantities. The second reason 
is that the equation of state more than any other 
quantity is sensitive to high momentum modes 
and therefore to cut-off effects. Indeed calcula- 
tions with standard staggered action as well as 
Wilson (clover) action show large dependence on 
the lattice spacing (Nf) '42" 33 . Using improved 
staggered fermions cut-off effects in thermody- 
namic quantities can be greatly reduced In 
Ref. |45j pressure and energy density have been 
calculated with p4 action for Nt — 4: and dif- 
ferent number of flavors and for quark masses 
corresponding to pseudo-scalar (pion) mass of 
about 700Af eV. When expressed in units of 
both pressure and energy density turn out to be 
much larger than in SU(3) gauge theory reflect- 
ing the presence of additional degrees of freedom 
(quarks). An interesting question is what is the 
energy density at the pseudo-critical point. Cal- 
culations with both Wilson and p4 staggered ac- 
tions yield consistent estimates, ec = (6 ± 2)T^ 
gSI, which is in turn not too sensitive to the 
quark mass and flavor. Assuming that this holds 
down to physical quark masses and taking the 
most recent value Tc = 169MeV one can esti- 
mate the (pseudo)critical energy density to be 
(0.6±0.3)GeV"//m3. 

In Fig. 01 we show the pressure normalized by 
continuum Stefan-Boltzmann limit for different 
number of quark flavors. One can see almost no 
flavor dependence in p/psB for T > Tc, i.e. the 
flavor dependence of the pressure in the high tem- 
perature phase can be well described by flavor de- 
pendence oipsB- Note that the pure gauge result 
was extrapolated to the continuum limit. This 
means that if the approximate flavor indepen- 
dence is true the p4 action results at Nt — 4 are 
already close to the continuum limit. Nonetheless 
one should keep in mind that the above calcula- 
tions were done at fixed niquark/T and not along 
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Figure 4. The pressure normalized by Stefan- 
Boltzmann limit for different number of flavors. 

the lines of constant physics. It has been shown in 
Ref. |47j that keeping niquark/T fixed instead of 
keeping the quark mass iriquark constant in phys- 
ical units can reduce the pressure by 9% at 2.5Tc. 

Recently the pressure has also been calcu- 
lated at non-zero chemical potential using both 
the Taylor expansion technique (38^ and multi- 
parameter re- weighting |47| . More precisely one 
calculates the baryonic contribution to the pres- 
sure Ap{T, = p(T, y) - p{T, A* = 0). The ratio 
Ap/ApsB, where ApsB is the ideal gas limit of 
Ap, is almost independent of the chemical poten- 
tial 38 47 . This means that the dependence of 
the pressure on the chemical potential is well de- 
scribed by that of ideal quark-gluon plasma. Tay- 
lor expansion method with p4 action gives similar 
result [2Sj which will be discussed in the next sec- 
tion. Here we only note that despite the similar 
behavior of Ap/ ApsB in the two methods there 
is a problem of cut-off dependence, in. Ref 02j 
ApsB is calculated for Nt = 4: lattice , while in 
Ref PHI the continuum value for ApsB was used. 

5. Comparison with the resonance gas 
model 

The QCD transition has two aspects : decon- 
finement and chiral symmetry restoration. In ab- 
sence of quarks in the fundamental representation 
an order parameter (the expectation value of the 
Polyakov loop) and a symmetry {Z{N)) can be 
linked to the deconfinement transition. In the 
more general case it is more meaningful to define 
the deconfinement transition as rapid increase of 
degrees of freedom. The chiral symmetry restora- 
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tion is linked to the behavior of the chiral conden- 
sate and Hghtest hadrons with increasing temper- 
ature. The interconnection between chiral and 
deconfinement aspects of the QCD transition is 
subject to intensive discussion in the literature 
|48|. There are at least two question which can 
be addressed in this respect. Does the mass of 
the lightest excitation (e.g. pion) determines the 
value of the transition temperature? What drives 
the QCD transition? The answer to the first ques- 
tion is probably "no" and can be explained using 
several reasonings. First, if the mass of the light- 
est excitation determines the value of the tran- 
sition temperature it would be difficult to un- 
derstand why the transition temperature in the 
SU(3) gauge theory, where the lowest excitation 
is the scalar glueball of mass of about 1.5GeV is 
only 50% larger than the transition temperature 
in full QCD where pions of mass lAOMeV are 
the lowest excitations. Second in terms of func- 
tional integral the finite temperature system can 
be viewed as a zero temperature system in a slab 
of extension 1/T (T being the temperature). The 
non-perturbative features of the vacuum (causing 
confinement and chiral symmetry breaking) are 
obviously lost when the size of the slab becomes 
too small. There is no obvious connection be- 
tween the "critical" size of the slab where the rel- 
evant non-perturbative features are lost and the 
mass of the lowest excitation ^ . 

In the following we are going to discuss the res- 
onance gas model which provides an answer to 
both questions raised above in terms of physi- 
cal degrees of freedom. Despite being a model 
it can explain many features of the available lat- 
tice data and can provide a standard reference for 
bulk thermodynamic properties in the low tem- 
perature phase similarly to the ideal quark-gluon 
plasma case ( Stefan-Boltzmann limit) which 
serves as a reference in the high temperature 
phase. 

In the resonance gas model the partition func- 
tion is a sum of one- particle partition functions 
of all mesons, baryons and their resonances |49j . 
In Fig. I show the lattice data for the en- 
ergy density (obtained with pA action) in com- 




Figure 5. The energy density in resonance gas 
model (thick line) against lattice data with p4 
improved staggered fermions j49| . 

parison with the resonance gas model which in- 
cludes about 2000 degrees of freedom correspond- 
ing to resonances with masses below 2GeV. Since 
lattice calculations are done at unphysical quark 
masses the masses of experimentally know reso- 
nances have been extrapolated upward in quark 
mass to take this into account 0^1 ■ Note that pi- 
ons contribute only 14% to the energy density at 
the transition temperature 49 . As one can see 
the resonance gas model can described the lattice 
data very well till temperatures of about l.lTc. 

One immediate consequence of the resonance 
gas model is the factorization of the baryonic 
pressure 

^ - P(^'^^-f^'^-'^ . FiT) cosh(|), (1) 

which holds as long {rrinucieon ~ 3/i) > T 50 . On 
the lattice one can calculate the pressure at finite 
chemical potential in form of the Taylor expan- 
sion 



^ C2{T){^r+c,{T)i^f + c,{T)i^ 



•(2) 



thank J. Kuti for bringing up this point. 



The factorization formula has important conse- 
quences for the baryonic pressure and its deriva- 
tives. The ratio of the baryonic pressure and 
its derivatives with respect to fj, (quark number 
and quark number susceptibility) are function of 
fi/T only The ratios of the expansion coef- 
ficients in Eq. jSJ are those of cosh: C4/C2 = 3/4, 
Ce/c4 — 0.3. These predictions can be compared 
with lattice data as shown in Fig. |S1 The 
figure clearly shows that the resonance gas model 
works quite well at low temperatures. 



7 




T/T„ 



Figure 6. The ratios of the expansion coefficients 
Cn calculated in resonance gas model 50 and on 
the lattice |5H . 
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Figure 7. The dependence of the deconfinement 
temperature on the pseudo-scalar mass calculated 
from lattice simulations [TSj and resonance gas 
model gni- 

Now we can return to the question what drives 
the transition in QCD. The energy density at 
the pseudo-critical point can be estimated to be 
Sc — 0.6GeV/ fm^ (see previous section). Assum- 
ing that this critical energy density is indepen- 
dent of the quark mass and number of flavors, 
as suggested by the existing lattice data, one can 
study the mass and flavor dependence of the tran- 
sition temperature Tc- This is shown in Fig. |3 
which shows that the resonance gas at fixed en- 
ergy density of O.SGeV/ fm^ seems to explain the 
quark (pion) mass dependence of Tc as long as 
iTips/y/o' < 4. For larger pseudo-scalar masses 
the discrepancies can be explained by the fact 
that at such large quark mass glueball states be- 
come important. It has been shown that inclusion 
of glueball states with thermal mass reduction can 
account for this discrepancy 

Thus we see that deconfinement happens once 
the energy density reaches a value of about 
O.SGeV/ fm^ and does not depend strongly on the 



details of the spectrum but only on its exponen- 
tial rise at large resonance masses. The fact that 
space gets (over)populated by hadron resonances 
near Tc also means that there is little space left 
for the non-perturbative vacuum if one assumes a 
simple MIT-bag motivated picture where the vac- 
uum inside a hadron is perturbative. This may 
lead to the decrease of the average value of chi- 
ral condensate. Also production of resonances in- 
creases the entropy, and the Polyakov loop, which 
can be related to the additional entropy generated 
by a static meson (see next section, for entropy 
increase due to adding an extra baryon into the 
system see Q2i ) , therefore rapidly increases near 
Tc. 

6. The free energy of static quarks 

It is customary to study in-medium modifica- 
tion of inter-quark forces, e.g. screening in terms 
of the free energy of static quarks separated by 
some distance r. Following McLerran and Svetit- 
sky the partition function of the system with 
static quark anti-quark pair (QQ) normalized by 
the partition function of the system without static 
charges can be expressed as correlation function 
of temporal Wilson lines 



Z{T) 

W{x) ^Vexp{- 1 
Jo 



(W{f}W\0)), 



1/T 



dTAo{x, t)). 



(3) 



(4) 



TrW is usually referred to as Polyakov loop. Tak- 
ing the logarithm of the above expression gives 
the free energy difference of the system with static 
quarks and without them at fixed temperature 
T. Since the quark and anti-quark can be either 
in color singlet or color octet state one should 
separate those contributions using the appropri- 
ate projection operators which gives the following 
expressions for the singlet and octet free energies 

exp(-Fi(r,T) + C) = ^{TiW {r)WHO)) , (5) 
exp(-F8 (r,T) + C) = -i-(TrT4^(f)T4^t(o)) 



4(TrTy(f)Triyt(o)). 



24 



(6) 
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One can define also the color average free energy 
which is the thermal averaged of the free energies 
in color singlet and octet channels |53I55| 

exp(-F(r, T) + C)^ UtvW [f)TvW\Q)) 
9 

= i exp(-Fi(r, T)) + ^ cxp(-F8(r, T)). (7) 

The normalization constant C in the above equa- 
tions will be defined below. The color averaged 
free energy is explicitly gauge invariant, while the 
others need gauge fixing or the Wilson line to be 
replaced by the gauge invariant Wilson line as this 
was done in Ref . ,56, . The results to be presented 
here were obtained using the Coulomb gauge as 
in this gauge a transfer matrix can be defined and 
the free energies have a meaningful zero tempera- 
ture limit. In fact it was shown that the Coulomb 
gauge and gauge invariant definition of the singlet 
free energy give numerically the same result ji6 . 

In Fig. |S1 we show the singlet, octet and 
averaged free energies in SU{i) gauge theory 
[57158159] . The normalization constant C was 
chosen such that the singlet free energies coin- 
cides with the zero temperature potential at the 
shortest distance. Below the deconfinement tem- 
perature all free energies are linearly rising, indi- 
cating confinement. The string tension, however, 
is smaller than in the vacuum. Above deconfine- 
ment all the free energies approach the same con- 
stant value at asymptotic distances indicating the 
screening in plasma phase. Thus the free energies 
of static QQ pair provide a true test of confine- 
ment. If this free energy is finite at infinite sep- 
aration the system is deconfined and otherwise 
it is confined. Usually the expectation value of 
the Polyakov loop (TrPF) is considered to be an 
order parameter for deconfinement and it is said 
that it is related to the free energy of a static 
quark \{TyW)\ ~ exp(-i^Q(T)/T). This quan- 
tity, however, does not have a meaningful contin- 
uum limit; for any temperature it vanishes in the 
continuum limit. One can define a renormalized 
Polyakov loop as Lren{T) = exp(-Foo(r)/(2T)) 
with Fao{T) being the asymptotic value of the 
QQ free energy at infinite separation (^T]. This 
provides an order parameter which is cutoff inde- 
pendent and has a well defined continuum limit 




Figure 8. The color singlet, octet and averaged 
free energy in SU(3) gauge theory above and be- 
low deconfinement 



|57j . While in the deconfined phase it is not dif- 
ficult to accept the concept of the octet free en- 
ergy, what it means in the confined phase is not 
quite clear. The construction in Eq. (j^ fixes only 
the relative color orientation of static charges, e.g. 
{QQ)s being in octet state. In order for such ob- 
ject to exist in the confined phase the net color 
charge must be compensated by gluons. Such an 
object can be identified with excited state of the 
string or hybrid potential |60I61| . In Ref. the 
definitions of the singlet and octet free energies 
were applied in the zero temperature limit. It was 
found there that the definition of the octet free 
energy does not project onto the string excita- 
tion as expected. This is because the amplitudes 
c„(r) in the spectral decomposition of the correla- 
tors J2n Cn('') exp(— i?„(r)/r) have non-trivial r- 
dependence At short distances where the r- 
dependence of c„ is likely to vanish and the octet 
free energy in © could give the hybrid potential. 

Notice that at short distances the singlet free 
energy coincides with the zero temperature po- 
tential in a certain distance range and is tem- 
perature independent, refiecting the fact that 
medium effects are not important at short dis- 
tances. This fact motivated our choice of the 
renormalization constant. To study the onset of 
medium effects and characterize the strength of 
interaction in the plasma it is convenient to in- 
troduce the effective running coupling constant 
Q!e//(r, r) = ''^dr''^^ which is shown in 
Fig. O At short distances the scale which de- 
termines the running of the coupling constant is 
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Figure 9. The effective running coupling constant 
for different temperatures j63| . 

given by the distance r and running follows the 
3-loop beta function. At larger distances the fast 
running of the coupling constant mimics the pres- 
ence of the linear term in Cornell parameteriza- 
tion of the zero temperature potential and there- 
fore is of non-perturbative origin. It is interesting 
to note that this non-perturbative feature of the 
running coupling constant survives also in the de- 
confined phase as shown in Fig. |51 The effective 
running coupling starts to deviate from the zero 
temperature result at some distance and eventu- 
ally decreases indicating the onset of medium ef- 
fects and screening. The distance where it hap- 
pens is obviously temperature dependent. Note 
that at distances where screening sets in, the sin- 
glet free energy does not have a simple screened 
Coulomb form expected in finite temperature per- 
turbation theory. Only at somewhat larger dis- 
tances the singlet free energy can be fitted by the 
form |as(T) exp(— mi(T)r)/r motivated by per- 
turbation theory • This defines then a cou- 
pling constant which depends on the temperature 
and not on the distance. In other words the run- 
ning of the coupling at finite temperature is con- 
trolled by the distance r for small r and by the 
temperature for r ^ 1/T. The results presented 
so far refer to SU (3) gauge theory. Very recently 
the study of the static QQ free energies was ex- 
tended to two and three flavor QCD 64 65 66 . 
The results for the singlet free energy in 3 flavor 
QCD are shown in Fig. ^| As in the case of pure 
gauge theory at short distances the free energy is 
temperature independent and coincides with the 




Figure 10. The singlet free energy in three flavor 
QCD "HS*. The temperature is given in MeV. 

zero temperature potential. However, the free en- 
ergy reaches a constant value also at low temper- 
atures, it flattens at distances of about O.Qfm. 
This is usually interpreted as string breaking. 

The free energies of static quark anti-quark 
pairs are often called and interpreted as poten- 
tials. They have been used in potential models 
to predict the pattern on J/ip suppression j67| . 
However, one should be aware that in fact one 
deals with free energies which contain an entropy 
contribution. Having the partition function in the 
presence of static QQ pairs one can also calculate 
the entropy and the internal energy difference be- 
tween the system with and without static charges 



Vi{r,T) = -T 



^dFi{r,T)/T 



Si{r,T) 



-T 



dT 

dFijr, T) 
dT ' 



(8) 



(9) 



The free energy is decreasing with temperature 
(see e.g. Fig. I1U|I . This implies that the entropy 
contribution is negative (c.f. equation above). In 
Fig. El we show the internal energy at differ- 
ent temperatures. It is obvious that the internal 
energy is larger than the free energy (compare 
with Fig. IHI ) . If one would interpret the internal 
energy as potential it would lead to the conclu- 
sion that the J /ip remains bound till 1.7Tc |68I69| . 
One may wonder to which extent such interpre- 
tation is correct. In Fig. ^Jwe show the asymp- 
totic value of the internal energy as function of 
the temperature. One can see a large increase 
of Uoo (T) in the vicinity of the transition. Such 
large increase is difficult to interpret in terms of 
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Figure 11. The singlet internal energy above de- 
confinement in SU(3) gauge theory 
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Figure 12. The asymptotic value of the internal 
energy energy in SU(3) gauge theory j^, 2 flavor 
QCD and 3 flavor QCD The fine is the 
asymptotic value of the free energy in 2 flavor 
QCD. 

modification of inter-quark forces and it is prob- 
ably due to many body effects. 

7. Meson spectral functions 

Lattice QCD can provide information about 
imaginary time correlation functions. On the 
other hand to make contact with observable quan- 
tities (e.g. masses and widths of hadron reso- 
nances in the medium) one should deal with cor- 
relators in real time 

D<{t) = (6(o)6(t)>, D>{t) = (6(t)o(o)>, (10) 

or equivalently the spectral function 
D>{lu)~ D<{lo) 



2-K 
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For example the rate of emission of thermal dilcp- 
tons can be expressed in terms of the spectral 



function of the correlator of the vector current 



dW 



5a 



1 



-av{po,p)- 



(12) 



Here we consider correlators of point meson op- 
erators qTq, with F = 1, 75, 7^,7^75 for scalar, 
pseudo-scalar, vector and axial-vector channels. 
We will also restrict the discussion to the case 
of zero spatial momentum. The imaginary time 
correlator is an analytic continuation of the real 
time correlator G(r) — D^{—iT). Using this and 
periodicity in the imaginary time (also known 
as Kubo-Martin-Schwinger condition) an integral 
representation for G{t) in terms of the spectral 
function can be derived 



G(r,T) 
K{t,uj) 



duJ<7{LU, T)K{t, lu, T) 

cosh{uj{T- 1/2T)) 
sinh(w/2T) ' 



(13) 
(14) 



This relation is valid in the continuum. It is not 
obvious that it holds also for the correlators calcu- 
lated on the lattice. However, it was shown that 
an integral representation like in Eq. 1141 can be 
also defined for the lattice correlator in the limit 
of the free field theory jTJ. The only modifica- 
tion is that the upper integration limit becomes 
finite Wmaa; — 4a~^ 71 . The cutoff dependence of 
the correlators is contained in the spectral func- 
tion ^L|. The obvious problem is that we have 
to reconstruct several hundred degrees of freedom 
(needed for a reasonable discretization of the in- 
tegral in Eq. ^^)- with 0(10) data point on 
G(t). This ill-posed problem can be solved using 
the Maximum Entropy Method (MEM) f?^. The 
method has been applied to study meson spectral 
functions at zero |72l73l74j as well as finite tem- 
perature IZSEaZZEHEMHU. 

The studies of meson spectral functions at zero 
temperature reproduced the properties of ground 
state mesons but on the other hand revealed 
structures (lattice artifacts) whose position scales 
like |73I81I74) . The cutoff effects in the me- 
son spectral functions extracted from the inter- 
acting theory have little to do with those in the 
free theory. 
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Figure 13. The charmonia spectral functions below deconfinement for two different lattice spacings 
a-i = 4.86Gey (left) and a"! = 9.72GeV (right) (H!]. 

tion may be quite small < 1 /T (see discussion in 
Ref. [7^1 on this point). While the former limita- 
tion can be relaxed using larger anisotropics the 
later will always remain. Therefore it is useful to 
inspect the temperature dependence of the corre- 
lators at the first step. The temperature depen- 
dence of the correlators is due to the explicit tem- 
perature dependence of the corresponding spec- 
tral functions and the temperature dependence of 
the integration kernel in Eq. (|14(l . To factor out 
the trivial temperature dependence of the corre- 
lator due to periodicity in the time direction we 
construct the following model correlator 



First I will discuss finite temperature meson 
spectral functions for heavy quarks. This is of 
particular interest as it was suggested by Mat- 
sui and Satz that existence or non-existence of 
heavy quarkonia in the plasma reflects the screen- 
ing properties of the plasma and can be used as 
signal for deconfinement and quark gluon plasma 
formation [H2|. Different calculations based on 
screening and potential models have estimated 
that J/il) will dissolve around l.lTc |67I83| . Cal- 
culations of quakonium spectral function using 
MEM were performed only recently. Calcula- 
tions have been done using anisotropic as well as 
fine isotropic lattices in quenched approximation 
|79l8UI8l |. Whatever approach is used several 
spatial lattice spacings are needed to check the 
cutoff dependence of the result. To demonstrate 
this in Fig. El we show charmonium spectral 
functions below deconfinement at two different 
lattice spacings from Ref. 81 . The pseudo-scalar 
and vector channels correspond to r/c and J /ip 
states, while the scalar and axial-vector channels 
correspond to Xco and Xci states. One can see 
that while the ground state peaks reproduce the 
expected masses the positions of the other peaks 
scales roughly as 1/a, i.e. these peaks are lattice 
artifacts. Thus point operators can provide infor- 
mation only about the properties of the ground 
state for a given channel. 

Now let me discuss what happens at finite tem- 
perature. In this case the reconstruction of the 
spectral function becomes more difficult as not 
only the number of available time-slices is reduced 
but also the extent of the imaginary time direc- 



dLoa{uj,T*)K{T,LU,T), 



(15) 



with T* being some temperature below deconfine- 
ment. If the spectral function has no temperature 
dependence also above deconfinement we would 
expect no temperature dependence in G/Grec and 
G/Grec — 1- In Fig. iniwe show G/Grec for dif- 
ferent channels. As one can see for temperatures 
T < 2Tc this ratio is not very far from unity for 
the pseudo-scalar and vector channels, while for 
the scalar and axial- vector channels large changes 
are visible already at l.lTc. This suggests that 
ground state charmonia may survive deconfine- 
ment with little modifications of their properties, 
while excited charmonia states are dissolved or 
strongly modified by the medium right after de- 
confinement. To study the problem more in detail 
one needs to reconstruct the spectral functions 
which are shown in Fig. 1151 The figure confirms 
the expectation based on the analysis of the cor- 
relator, J/ip seems to exist up to temperatures 
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Figure 14. Ratios G/Grecon for the pseudo-scalar and the vector correlators (left) and for the scalar and 
the axial- vector correlators (right) [SJj . For better visibility in the pseudo-scalar and axial- vector channel 
the data for G/Grecon have been shifted by a constant. 
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Figure 15. Charmonia spectral functions for the vector (left) and scalar (right) channels above Tc- |81) . 



2.25rc , while the Xco states is dissolved at l.lTc. 
The situation is similar for the pseudo-scalar and 
axial vector spectral functions. The results for 
the IS' states {J/ip, rjc) are confirmed by calcu- 
lations done on anisotropic lattices, though there 
is some controversy about the temperature where 
J/ip dissolves and the way it disappears. In Ref. 
|80| it was found that the IS" charmonia abruptly 
disappear at l.TTc, while the spectral functions 
shown above suggest a gradual dissolution. 

Meson spectral functions at finite temperature 
were also studied in the case of light quarks. In 
particular in Ref. [75] the first attempt to cal- 
culate vector spectral functions and the thermal 
dilepton rate was made using isotropic lattices. 
Then meson spectral functions for scalar, pseudo- 
scalar, vector and axial- vector channels were cal- 
culated using anisotropic lattices (SO] . One should 
keep in mind, however, that unlike in heavy quark 
case here no study of possible systematic effects 
was done so far. One common expectation that 
the pseudo-scalar and scalar correlators should 
become degenerate at high temperatures due to 
effective restoration of C/4(l) symmetry. The 



same is true for the vector and axial- vector corre- 
lators. Explicit calculations confirm these expec- 
tations |76l78j . An interesting feature of the light 
meson spectral function above deconfinement is 
that they do not show any continuum like struc- 
tures which one could identify with freely propa- 
gating quarks but peaks. In Fig. ^|I show the 
meson spectral functions at lATc which demon- 
strates this feature. While the second and third 
peak are likely to be lattice artifacts (though not 
yet shown by a detailed analysis) the first peak 
is likely to be physical and its existence is quite 
bizarre. Note that the spectral functions in Fig. 
[TBI also show the degeneracy of opposite parity 
channels as expected from the analysis of the cor- 
relators, moreover the masses in all four channels 
appear to be roughly the same. 
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